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We perform experiments on thin-film diaphragms to show that the folding 
patterns of anisotropically compressed diaphragms are strikingly different from 
those of isotropically compressed ones. We then use a simple von Karman 
model to relate the overall features of these folding patterns to the underlying 
energetics. We show that the differences between the isotropic and anisotropic 
cases can be traced back to fundamental changes in the energy structure of 
the diaphragms. Finally, we point out that the energy structure of thin-film 
diaphragms is similar to that of many other systems in physics and engineering, 
into which our study may provide interesting insights. 

Consider a film of thickness h bonded to the flat surface of a substrate except over a 



domain Q of characteristic size d, Fig. |](a). The portion f2 of the film is a diaphragm 
[1] . If we apply a strain e* to the substrate, in the plane of the film, the diaphragm may 
deflect out of the substrate and fold [2] . We study the conditions which e* must fulfill for 
this to occur, and the nature of the resulting folding. 

We prepared diaphragms by gluing paper sheets and polymeric films (h ~ 0.01 to 
0.1mm) onto substrates (thickness ~ 30 mm) of the shape shown in Fig. |l](a), made of a 
high-density styrofoam [3]. Then we compressed the substrates in two perpendicular direc- 
tions, x\ and X2, using screw-driven steel plates. When the applied strain components are 
equal to each other, e\ = e%, the strain is isotropic; otherwise it is anisotropic. Figs. |](b-c) 
show the folding of two equally shaped diaphragms subject to isotropic strains. The fold- 
ing pattern is the same in both diaphragms; the number and spatial arrangement of the 
folds depend exclusively on the shape of being independent of d/h and of the strain. 
A completely different situation obtains when the strains are anisotropic. The folds are 
then perpendicular to the direction of e\ (where e\ > e\ by convention), regardless of the 
shape of Q, Figs. |T](d-e). The number of folds depends on d/h, however, see Figs. [l|(e-g); 
a comparison of Figs. |l|(e) and (h) reveals that it depends also on the strain. 

For the analysis of these results we follow [2] in i ) modeling the film as a von Karman 
plate, and ii) constraining the in-plane displacements of the diaphragm to remain null. 
As we shall show, this simple model allows for a straighforward energetic interpretation of 
the overall characteristics of folding patterns in thin-film diaphragms; the study of specific 
features of the folding, as well as of folding in structures other than diaphragms, may 
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require the consideration of the in-plane displacements [4]. For the proposed model the 
bending energy density at a point , x,2 ) of the diaphragm is 
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and the membrane energy density is 



V =2 



1 {\V W \ 2 -2{e{+ve* 2 )f + {l-v)w 2 (ej 



+ <in, (2) 



_4 

where C = Eh/(1 — u 2 ) is the membrane stiffness of the film; E is the Young modulus 
and v the Poisson ratio of the film; w is the out-of-plane deflection of the diaphragm; 
Vto = (w,i ,w,2 ) is the gradient of w; \Vw\ = (w 2 +w, 2 ) 1 ^ 2 is the largest slope of the 
diaphragm at the given point; and we have defined V'min = ~ ^ 2 )( e 2) 2 /2- The total 
energy of the diaphragm is the sum of the bending and membrane energies, $ = <3? m + $ b , 
where $ m and <3? b are the integrals over Q of ip m and ip h , respectively. We search a folding, 
described by w(xi, X2), which minimizes <&[w] subject to the conditions w = and w, n = 
on the boundary, Fig. |l](a). 

We start by noting that $ b and $ m are of order h 3 and h, respectively. Since we are 
interested in the thin-film limit, h — > 0, we expect <3? m to admit one or more minimizers 
with the same overall structure as the minimizer of We know that in the minimizers 
of <I> m the folds will not be rounded, but take the form of lines of slope discontinuity or 
sharp folds. (The reason is that ( £ m contains first derivatives of w only.) The bending 
energy will be confined to these sharp folds, in the form of an energy per unit length of 
sharp fold [5,6]. We plan to i) find the foldings which minimize $ m , and then ii) select 



among them the one which contains the least sharp-fold energy; this one we shall call the 
preferred folding [2,7,8]. 

We try to minimize <3? m by minimizing its integrand, ip m , Eq. |2[ For convenience 
we define a compressive regime e\ + ve\ > 0, and a characteristic slope k = ^Je\ + ve\. 
We look for gradients Vro = (u>,i , w,% ) for which (p m attains a minimum. When e* 
falls outside the compressive regime, a single minimum exists for Vtw = (0,0), and the 
diaphragm remains flat. In the isotropic case, = e^, infinitely many minima of value 
^min ex ist i n t ne compressive regime; they occur for |Vw;| = k, i.e. whenever the largest 
slope equals the characteristic slope. In the anisotropic case, two minima of value </?™ in 
exist in the compressive regime; they occur when Vw = (±/c, 0). Thus for a diaphragm in 
the compressive regime the minimum possible of the membrane energy is <]?™ f = 9?™ in ^4Q, 
where Aa is the area of the diaphragm. For convenience we shall work with y? m = (p m — 
<^™ in , which implies |> m = $ m - $PJ f) and l>P f = 0. 

We study the isotropic case first. Consider a ID example in which we impose the 
constraint w,2 = 0. Then, the minima of </3 m occur for w,i = ±k, and we can construct 
a minimizer of $ m by covering £1 with any set of simple roofs of slopes ±k (such as R\, 
i?2 and R3 in Fig. ||(a)), and then choosing their upper envelope. Infinitely many such 
minimizers exist, all of which contain sharp folds, as expected. To minimize the energy 
associated with the sharp folds, we must minimize their number; this we effect by selecting 
the upper envelope of all the minimizers, Fig. ||(a). In 2D a preferred folding can be found 
analogously, as the upper envelope of all the possible cobertures of £1 with conical roofs of 
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slope k, Fig. |2|(b) [2,9]. The preferred folding so selected does i) minimize <3? m ; ii) depend 
exclusively on the shape of f2; and Hi) match the experimental observations, as shown by 
a comparison of Fig. ^(c) with Figs. 0(b) and (c). 

Now we turn to the anisotropic case. Consider first the infinite diaphragm of Fig. ||(a). 
For this diaphragm the preferred folding can be found as in the ID isotropic case, in 
the form of a single simple roof of slopes w,\ = ±k, Fig. ||(a). In the finite diaphragm 
of Fig. ||(b), on the other hand, a single simple roof violates the boundary conditions 
w{x2 = 0) = w(x2 = 3/2) = 0. To circumvent this problem we try the folding wq, of 
Fig. ||(b). On the triangular regions of Fig. ||(b), which we call closure domains, w,\ ^ ±fe, 
and therefore y>™ > fi^m- It follows that wo is not a minimizer. In fact, calling ^4a[^o] the 
area of the closure domains, $ m [u;o] = ^a^a^o] > ^Sf We can a PP roacn the minimum 
value <!>™ f = 0, however, by using foldings Wj with j = 1,2,..., e.g. Fig. ^C, for which 
^A^j] ~^ and therefore 3> m [w;j] = (/j^^a^'] — * as j — > oo. Thus <£> m can be made 
arbitrarily close to its minimum by allowing the diaphragm to become highly folded. The 
sequence Wj is said to be a minimizing sequence, and the associated foldings are called 
micro structures [10,11]. 

For arbitrarily shaped diaphragms, we can construct minimizing sequences in an anal- 
ogous way, e.g. Fig. ^(d). All the microstructures in a minimizing sequence contain sharp 
folds and closure domains. As the microstructures become more folded, the sharp-fold 
energy increases, and the closure-domain energy decreases. In principle, we can identify a 
preferred microstructure for which the trade off between closure-domain energy and sharp- 
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fold energy is resolved in the least total energy. Thus bending checks the infinite folding 
implied by a minimizing sequence [10]. A straightforward dimensional analysis reveals 
that the number of folds N in the preferred microstructure scales with ip^d/T, where T 
is the energy per unit length of sharp fold. Setting e\ = for the sake of simplicity, it is 
easy to conclude that </?™ oc E(e\) 2 h\ a more involved analysis leads to T oc E(el) 3 / 2 h 2 
[12]. Therefore, 

iVoc^)^. (3 ) 

According to this expression, the number of folds for the diaphragms of Figs. |l](d-h) should 
be proportional to 32, 30, 16, 10 and 22, respectively. The observed numbers are 32, 28, 
24, 12, and 22 (folds spanning the whole diaphragm). We conclude that the preferred 
microstructure of the minimizing sequence does i) fail to minimize $ m ; and ii) contain, 
in accord with the experimental evidence, folds which are perpendicular to the direction 
of e*, regardless of the shape of f2, in a number which scales with d/h and e* in the form 
ofEq. (|). 

To sum up, the folding of compressed diaphragms can be interpreted in terms of 
two operations: fold to release membrane energy, and then allow bending to select one 
among many possible foldings. In the isotropic case, infinitely many foldings exist which 
can accomodate the boundary conditions and simultaneously minimize & m . Out of these 
foldings bending selects a preferred one. In the anisotropic case, no folding exists which 
can accomodate the boundary conditions and simultaneously minimize & m . It is possible, 
however, to construct sequences of increasingly fine foldings, or microstructures, whose 
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associated membrane energies converge to the minimum value of $ m . Out of these foldings 
bending selects a preferred one. 

The occurence of microstructures in compressed diaphragms discloses a mathematical 
similarity with many problems including solidification [13], solid-state phase transforma- 
tion [14], epitaxial thin-film growth [15], crystal plasticity [16], ferromagnetism [17], et 
cetera, where microstructures have been documented in the form of eutectic structures, 
twinning, film roughening, dislocation cells, and magnetic domain structures, respectively. 
Other such problems are coagulation [18], stretching of solid foams [19], and self-assembly 
of polymer layers on patterned substrates [20]. Because experiments can be easily per- 
formed on compressed diaphragms, and a simple model appears to explain their behavior 
well, the study of diaphragms may prove useful to gain insights into many other systems. 
For example, a 90 degree rotation of the gradients in Figs. |3|(b-d) leads to closed-flux 
vector fields such as they obtain in ferromagnets, where magnetic poles are energetically 
penalized [17]. In this analogy, sharp folds model ferromagnetic domain walls, and strain 
anisotropy models crystalline anisotropy. A comparison of the fold branching observable 
close to diaphragm boundaries, e.g. Fig. 0(g), with the analogous phenomenon of domain 
branching in the vicinity of free surfaces in ferromagnets [21] would enhance our under- 
standing of the underlying energetics. Supported by a grant from the Mechanics and 
Structures of Materials Program, NSF, Dr. K. P. Chong, Program Director. 
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Figure 1: a, Thin film/substrate system, diaphragm, and notation. Top view of folded 
diaphragms: b, SiC film, e\ = e* 2 = 0.011, d/h = 140 [22]. c, Paper film, e\ = e* 2 = 0.040, 
d/h = 950. The following are polymeric films (except when noted) with e?i = 0: d, e\ = 
0.038, d/h = 5600. e, e{ = 0.035, d/h = 5600. f, e\ = 0.035, d/h = 2840. g, Paper film, 
el = 0.040 d/h = 1800. h, e\ = 0.019, d/h ^©600. 



Figure 2: Analysis of the isotropic case, a, Side view of a ID, constrained diaphragm 
covered with a set of three simple roofs of slopes ±/c, the associated minimizer (bold lines), 
and the preferred folding (dashed lines), b, A 2D diaphragm and an example of conical 
roof of slope k, denoted C\. c Preferred folding for the diaphragm shape of Figs. [I|(b) and 
(c). 




a b c d 



Figure 3: Analysis of the anisotropic case, a, Top view of a semi-infinite diaphragm with 
preferred folding, b, Rectangular diaphragm with folding Wq. c, Idem with u>2 (i.e. the 
microstructure j = 2 of the minimizing sequence Wj). d, Two terms of a minimizing 
sequence for a circular diaphragm. 
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